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Electronic dispersion anomalies in the iron 
pnictide superconductor Bai_ x K x Fe2As2 
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Institut fur Theoretische Festkorperphysik, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany 
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3 Fachbereich Physik, Universitat Konstanz, D-78457 Konstanz, Germany 

The pairing mechanism in the iron-pnictide superconductors is still unknown. How- 
ever, similarities to the cuprate high-temperature superconductors suggest that a simi- 
lar mechanism may be at work. Recently, careful experimental studies of the spin exci- 
tation spectrum revealed, like in the cuprates, a strong temperature dependence in the 
normal state and a resonance feature in the superconducting state. Motivated by these 
findings, we develop a model of electrons interacting with a temperature dependent 
magnetic excitation spectrum based on these experimental observations. We apply 
it to analyse angle resolved photoemission and tunnelling spectra in Bai a; K a: Fe2As2. 
We reproduce in quantitative agreement with experiment a renormalisation of the 
quasiparticle dispersion both in the normal and the superconducting state, and the de- 
pendence of the quasiparticle linewidth on binding energy. We estimate the strength 
of the coupling between electronic and spin excitations. Our findings support the 
possibility of a pairing mechanism based dominantly on such a coupling. 



Shortly after the discovery of high-temperature super- 
conductivity in Fe-based pnictide compounds [l|]a mag- 
netic Cooper-pairing mechanism was proposed [2j , while 
electron-phonon interaction as primary pairing mecha- 
nism was found to be unlikely jlj. This conjecture is 
supported by the proximity of antiferromagnetism and 
superconductivity in the phase diagram |4j-l8(. The spin 
excitation spectrum in pnictides shows pronounced sim- 
ilarities with other superconductors where a magnetic 
pairing mechanism is under debate. In particular, the 
strong temperature dependence of the normal state spin 
excitations studied recently by Inosov et al. 0], as well 
as the presence of a spin resonance feature in the super- 
conducting state [Io| are prominent features also present 
in cuprates, as well as in some heavy fermion supercon- 
ductors. 

Furthermore, angle resolved photoemission (ARPES) 
measurements reveal a sharp Fermi surface consisting of 
electron-like and hole-like pockets that exhibit compara- 
ble superconducting order parameter amplitudes and are 
nearly nested by the antifcrromagnctic wave vector Q 
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12f , a reason to believe that magnetic and electronic 

This raises the ques- 
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order are closely connected 
tion how strongly electrons couple to spin fluctuations, 
as those lead to an effective electron-electron attraction, 
if the order parameter changes the sign on the differ- 
ent pockets 14 - IB| . Investigating the low energy disper- 
sion anomalies, whose position and shape can be traced 
back to a coupling of bosonic modes, has proven to be 
a p owerful method for obtaining this information 
Il9j . Recent experimental studies on the superconducting 
state spectral function of Bai_ a; K 2 ;Fe2As2 measured by 
ARPES find such an anomaly, which appears as a kink 
in the dispersion relation at about 25 meV [l|| . On the 



other hand, inelastic neutron scattering (INS) studies of 
this compound reveal the development of a resonant spin 
excitation in the superconducting state [l(J. The tem- 
perature dependence of both neutron intensity and the 
self-energies extracted from ARPES follow an order pa- 
rameter like evolution. This strongly indicates that the 
spin resonance indeed accounts for the observed anomaly. 

Below, we present a theoretical model that explains 
the anomalous features in the dispersion relation by a 
coupling of electrons to a spin fluctuation spectrum that 
(a) reproduces the experimentally observed temperature 
dependent spectrum in the normal state and (b) features 
a low energy resonance in the superconducting state. We 
show that coupling to such a spin-fluctuation spectrum 
leads to a renormalisation of the quasiparticle dispersion 
even in the normal state and can account for the renor- 
malisation factor necessary to match density functional 
theory (DFT) calculations to experiment. Furthermore, 
taking into account the resonance in the superconducting 
state, we are able to quantitatively reproduce experimen- 
tal ARPES results, including low-energy and high-energy 
renormalisation of the dispersion and the energy depen- 
dence of the quasiparticle linewidth. Our theoretical SIN 
tunnelling spectra obtained for the same parameter set 
are consistent with recent experimental observations that 
show an asymmetric lineshape of the tunnelling conduc- 
tance [2p| . 

Spin fluctuation spectrum and resonance mode 

The appearance of a resonance in the dynamic magnetic 
susceptibility upon entering the superconducting state 
is well known from cuprate superconductivity. It is sit- 
uated at an energy below the particle-hole continuum 
and is peaked around the antifcrromagnctic wavevector. 
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FIG. 1: Fermi surface geometry and characteristics of the spin fluctuation mode, (a) The dominant orbital contri- 
butions to the bare Fermi surface at k z — —n, —n/2, 0. Red, green, blue and yellow correspond to the d xz , d yz , d xy , d x 2_ y 2 - 
orbitals, respectively, (b) Momentum dependence of the spin fluctuation mode. The mode is centred at wavevectors 
±Qi = (0, ±7r) and ±Q2 = (±7r, 0). (c) Energy-momentum dependence of the spin susceptibility in the normal state at 
T = 50 K (bottom) and in the superconducting state at T = 15 K (top), (d) The shape of the spectral function x"(uj,Q) 
is shown for T = 50 K (blue line) and T = 15 K (red line). The dashed green line shows the behaviour for T =15 K if the 
material were in the normal state, (e) The ratio of the resonance and the continuum spectral weight, wt/x.t as function of 
temperature. Concerning temperature dependence of x' ( w > 0) see Fig. [7] 



Although there arc various theoretical models for the 
spin-spin response function, here we prefer to not rely 
on a specific theoretical model, but rather use a semi- 
phenomcnological approach where the spin susceptibility 
is taken from experiment. This approach was very suc- 
cessful in the case of cuprates 17J , and has the advantage 
of being independent of specific theoretical assumptions 
and approximations. In addition, most of the magnetic 
correlations are included automatically within this ap- 
proach by relying on experimental results. 

A detailed investigation of the spin dynamics in the 
normal state and the temperature dependence of the spin 
resonance in pnictides was recently performed on opti- 
mally doped BaFei.85Coo.i5As2 by Inosov et al. @- They 
observe that the imaginary part of the normal state sus- 
ceptibility near the antiferromagnetic wave vectors Q a 
obeys in good approximation the Ornstein-Zernicke form 



-XT 



J? 2 l + # \q-Qa\' 



(1) 



where \t = Xo/(T + 9) controls the strength of magnetic 
correlations, ^ lax = £Iq(T + 9) sets the typical spin- 
fluctuation energy scale, £t = Co/ ^/T + 9 is the magnetic 
correlation length and 9 is the Curie- Weiss temperature. 



With decreasing temperature, spectral weight is 
shifted towards lower energies. Below T c the spectrum 
becomes gapped and a resonance appears at an energy 
that follows approximately an order parameter like evo- 
lution, i.e. fl^ es = \/l — T/T c . In order to model the 
resonance below the quasiparticle continuum we assume 
a (sharply peaked) Lorentzian in energy and a separable 
form for the momentum dependence, 



Xr,n(w,qr) 
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with the weight functions 



- S. r 



& \q - Qc 



(2) 



(3) 



In order to take care of the periodicity in reciprocal space 
we replace the factors \q — Q a \ 2 in Eqs. (JlJ and ((3|) 
by 4 [sin 2 (3kz^-) + sin 2 ( gi/ ~ 2 q ° !/ )], and we neglect the 
dependence on q z for our purpose as it varies weakly 
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22} . The resulting momentum dependence at the 
resonance energy is shown in Fig. QJb) , with the two an- 
tiferromagnetic wavevectors Qi = (n, 0) and Q2 = (0, tt). 
These wavevectors connect Fermi surface sheets with 
same orbital character (n = 1,2- • -5 corresponding to 
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TABLE I: Parameter Set used for Bao.6Ko.4Fe2As2 



o 



6) 



XT=OK 

15.6 meV 0.375 SgV 1.07 K 1/2 30 K 2 15.5 meV 3 meV 



the five Fe orbitals d xz , d yz , d xy , d x i_ y i, d 3z 2_ r 2), as 
illustrated in Fig. QJa) . In our model we consider only 
the main orbital contributions to quasiparticle scatter- 
ing resulting from the n = 1,2,3 orbitals and describe 
the different momentum variations with the help of the 
parameter b n>a (= 1/2 for na = 11,22,31,32 and = 
else). 

For the gapped continuum (in the absence of detailed 
experimental data) we make the approximation that for 
ui above a temperature dependent continuum threshold, 
ui > 2A(T), it is identical to its normal state value, 
Eq. ([1]), at the respective temperature. We note that 
this approximation is excellent except possibly close to 
the continuum threshold. Here, A is the mean supercon- 
ducting gap at the nested Fermi surfaces. Thus, denoting 
the imaginary part Im% of the retarded dynamical sus- 
ceptibility with x" i we approximate the superconducting 
state spectrum of the susceptibility by 



X"c.n(^<2) = 



X?,„(w,g) for \u\ <2A 
X " {u,q) for M >2A 



(4) 



and calculate the real part by exploiting Kramers-Kronig 
relations. In order to reduce the number of parameters 
we employ a local sum rule to determine the ratio be- 
tween resonance and continuum part, wt/xt- It is cho- 
sen so that the energy and momentum integrated spin 
structure factor, J™ du f dqS(q,u), with S(q,uj) = 



2ftx"(9,w)/(l-e 



-hu)/k B T 



) remains temperature inde- 



pendent in the normal as well as the superconducting 
phase (for technical details see the Methods section). The 
overall weight xo w ih be combined with the coupling con- 
stant g below, and it is only the combined quantity g 2 Xo 
that is of relevance for our results. However, as discussed 
below, from experimental values of xo an( 4 horn a de- 
termination of the fit parameter g 2 Xo an estimate of the 
coupling constant can be obtained. 

We apply this model to hole-doped Barj.6Ko.4Fe 2 As2, 
where data are available both for spin excitations as well 
as for electronic excitations, and where also the most 
detailed experimental angle resolved photoemission data 
exist. In FIG. Q] the momentum and energy dependence 
for the parameter set in TABLE Q] is presented. The 
resonance in the spin excitation appears at an energy 
&res 7K ~ 14meV below the quasiparticle continuum 
uj < 2A(15K) w 24meV [lfj. In momentum space the 
mode is peaked around the wave vectors Qi = (tt,0, q z ) 
and Q2 = (0,Tr,q z ) with correlation length £ r of nearly 
twice the lattice constant [FIG. H^b)]. The momentum 
dependence in z-direction is assumed to vary weakly as 



motivated by [21 , 22 1 . Comparing the energy-momentum 
distribution in the normal and the superconducting state 
[FIG. Hfc)], we see that low-energy spectral weight is 
shifted into the resonance when the continuum gap 
opens below 2A. In addition, the sharpened momentum 
distribution leads to an enhancement of the resonance 
spectral weight. In FIG. [He) we show the numerically 
calculated ratio wt/xt, which features this evolution 
and is in excellent agreement with the functional form 
wt/xt = w t =ok/xt=ok (1 - T 2 /T 2 ). As we will see 
later on, such a modelling will lead to an effect on the 
electronic dispersion that fits well with experimental 
observations. 

Coupling to spin fluctuations 

We are interested in the renormalisation of the fcrmionic 
dispersion as a result of the coupling of electrons to the 
spin fluctuation mode. The idea is to extract the influ- 
ence of the resonance by comparing the superconducting 
and normal state dispersion from which self-energy ef- 
fects are directly inferred. This allows for an immediate 
comparison with experiment. Following the approach for 
the cuprates 17J, one could assume that the momentum 
dependence in Eq. ((T|) and Eq. ((2j) exclusively chooses 
the Fermi surface sheets that are coupled by the mode. 
However, fluctuation exchange (FLEX) approaches have 
shown that the magnetic mode predominantly scatters 
between states with the same orbital character 27], |28( . 

To take this into account, we employ a tight-binding fit 
in orbital basis which was obtained from the DFT band 
structure of BaFe2As2 by Graser et al. (2f|. Details are 
discussed in the Methods section. 

In FIG. [T^a) we show the bare Fermi surface for k z = 
—it, — 7r/2 and 0, corresponding to a be tetragonal unit 
cell and the I4/mmm symmetry of the crystal. The hole 
pockets at (0, 0, k z ) and (71", 7r, k z ) are nearly nested to the 
electron pockets at (0,tt, k z ) and (7r,0, k z ) by the wave 
vectors Q\ and Q 2 , even when taking into account the or- 
bital characters [25l428| . We couple electrons to the spin 
fluctuation spectrum with a coupling constant <?, which 
is assumed to be independent of energy, momentum and 
orbital number. 

In ARPES experiments the intensity of photo elec- 
trons is proportional to f(e)A(e,k), where / is the 
Fermi distribution function. Dispersions obtained 
from those usually differ from DFT calculations by an 
renormalisation factor of 2 2^-3l|. In the upper panel 
of FIG. [2] we present the influence of a coupling to spin 
fluctuations in the normal state. It shows an intensity 
plot of the spectral function A(e, k) along a cut in the 1 st 
Brillouin zone, obtained numerically as explained in the 
Methods section. The band maxima and minima as well 
as the Fermi velocity very well align with experimentally 
observed values [2^]. By choosing the coupling constant 
g to be the same for all orbitals the band curvature 
becomes strongly enhanced leading to the observed 
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FIG. 2: Electronic spectral intensity. On top: spectral function in the normal state for T = 50 A" along a cut in the 1st 
Brillouin zone for k z — 0. At bottom: the same for the superconducting state at T — 15 K for the square regions indicated (a) 
and (b) in the top panel. The black lines show MDC-derived dispersions for the (a) u\ and a.2 band and the (b) f}\ band. At 
the lower right the renormalised Fermi surfaces including the group notations are shown. 



shallow electron pocket at the X-point. One clearly 
sees that the bands are bent by a factor 1.5 — 2 and 
that the quasiparticlc width increases with energy. Thus 
coupling to spin fluctuations gives an essential contri- 
bution to the above mentioned rcnormalisation factor. 
However the Fermi surface remains nearly unaffected as 
can be seen by comparing the Fermi surface at 50K on 
the right bottom of FIG.|2]with the bare one in FIG.QJa). 



Superconducting order 

The origin of the pairing instability may well be re- 
lated to the spin fluctuation continuum, as demonstrated 
by a recent FLEX calculation .23]. Our model is re- 
stricted to the low energy region in the spin excitation 
spectrum, and does not precisely treat the incoherent 
high-energy part. However, this part considerably con- 
tributes to pairing, whereas the energy range of interest 
here (|w| < u> c = 200 me V) gives only a partial contribu- 
tion (about 40-50%) to the value of the superconducting 
order parameter. Thus, we add in our theory a contri- 
bution Afc to the order parameter that results from the 
incoherent high-energy part of the spin-fluctuation spec- 
trum beyond 200 meV. 

The order parameter is chosen to have an s^-symmetry 
(here and in the following the unit of length is the in- 



plane lattice constant a) 



A fc (T) = A (T) cos{k x ) cos{k y ). 



(5) 



This pairing state is supported by experiment as well 
as numerical calculations [ll|, [HI The mag- 

nitude of the superconducting gap was observed to be 
A(15K) w 12meV at the inner hole-like pocket as well as 
the corresponding nested pockets HHI2. We choose A 
in Eq. ((SJ) so that the renormalised gap reaches the ex- 
perimentally observed value at these particular points in 
the Brillouin zone, which gives Ao(0) = 18.1 meV. The 
renormalised gap is obtained by taking into account a 
renormalisation factor as discussed in the Methods. For 
the temperature dependence (in the absence of detailed 
data) we assume the form A (T) = A (0)^/1 - T/T c . 

Dispersion and linewidth: self-energy effects 

In the superconducting state the dispersion features in 
FIG. [^a) are modified due to (a) the appearance of the 
superconducting gap, and (b) the modifications in the 
spin excitation spectrum that is coupled to the conduc- 
tion electrons. In the lower panel of FIG. [2] extracts from 
the hole-like and electron-like pockets at the V- and Ap- 
point are presented for 15 K, well in the superconducting 
state. At first glance the intensity plots show no clear 
hint of the bosonic resonance, in contrast to the eye- 
catching break features in the cuprates 32[. In order 
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to extract a quantitative effect we need to analyse dif- 
ferences between the dispersions and lincwidths in the 
normal and superconducting state in detail. 

Because in the case of pnictides multiple orbitals are 
involved in the electronic spectra, the linewidth function 
and the renormalisation of the dispersion cannot be re- 
lated in an easy way to theoretically obtained self energies 
as it has been possible in the case of the cuprates. For 
this reason we follow here a different path and extract 
quantities from theoretically obtained spectral functions 
exactly as they are extracted in the ARPES experiments 
(lil [3| from measured spectra. For fixed energy the 
momentum dependence of the spectral function (a so- 
called momentum distribution curve, or MDC, see inset 
to FIG. HJ is peaked, with the peak often well approxi- 
mated by a Lorentzian. If the self energy does not vary 
much as a function of momentum over the width of such 
an MDC peak, the MDC spectral function is of the form 



1 



A t (k) = - 



y" 



7r{« 6 (fc-fc e )} a + {E" y 



(6) 



where k e determines the maximum point of the MDC, 
and v e = {dK/de)- 1 0. In the bottom panel of FIG. M 
we show dispersions obtained from MDC maxima of the 
ai,2 and the /3i band as black curves. Furthermore, we 
obtain the linewidth in analogy to experiment by fitting 
a Lorentzian to the theoretically obtained ARPES 
spectra. 

We now discuss the thus obtained quantity S"(e) = 
that determines the linewidth of the MDC. Assum- 
ing that it weakly depends on momentum in a small 



region around fc c it is given by £ e ' 



v f 8k f , and 



|5fc e | is the full width at half maximum (FWHM) of the 
Lorentzian in reciprocal space in direction of v e . We un- 
derline that it is important to distinguish between the 
calculated self-energy S^(e) for the n th orbital, as dis- 
cussed in the Methods, and the experimentally motivated 
quantity S"(e) extracted by fitting a Lorentzian to the 
theoretically obtained ARPES spectra. Since the Green's 
function matrix in orbital space cannot be inverted an- 
alytically there is no simple correspondence between the 
two quantities above. We will call the experimentally 
motivated quantity S(e) imaginary part of the effective 
self energy below. 

As can be seen in FIG. [3] the imaginary part of the ef- 
fective self-energy exhibits a linear dependence for higher 
energies, i.e. S"(e) cx e, consistent with marginal Fermi 
liquid theory |34| . This results from the coupling to the 
continuum, in particular from the slow decay of the spin 
fluctuation spectrum towards high energies. We note 
that the correct magnitude of the linear in e high-energy 
part of £"(e) restricts the coupling strength g 2 Xo- We are 
able to reproduce the experimental observations [l9| with 
a coupling strength g 2 Xo = 1-17 x 10 3 /x|cV K. With this, 
all parameters of the theory are fixed by the high-energy 
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FIG. 3: Scattering rate: Imaginary part of the effective self- 
energy at T = 15 K determined by the width of the Lorentzian 
fits (centred at the points indicated in the inset) times the 
group velocity at the given energy, i.e. S" = v e 8k t , at the ai 
(red) and the Q2 (green) band. For temperature dependent 
results see Fig. [8] 



behaviour of the electronic spectra, and consequently the 
magnitude of all low-energy features are true predictions 
of the theory. 

The spin resonance leads to a hump feature in S"(e) 
for excitation energies between 20 — 30 meV at the ot\- 
band, which is apparently observed in experiment [l9| . 
This feature is washed out in the a2-band. However the 
actual position of the resonance can not be obtained from 
such a vague energy range. In order to give a convincing 
quantitative statement, one has to consider the modifica- 
tion of the dispersion upon entering the superconducting 
state, which we do now. 

The values k e determine an MDC-derived dispersion 
e = £k, as can be seen as black curves in FIG. Eta) and 
(b) . The modification of the bare dispersion to Ek can 
be expressed in terms of the real part of the effective self- 
energy, £', that is defined by £k = £,k + fc). The 
bare dispersion is not a measurable quantity. However, 
as the bare dispersion is temperature independent, any 
changes of £'(£&,&) = T,' k with temperature can be ex- 
tracted experimentally by taking the difference between 
two MDC-derived dispersions at fixed fc, provided again 
the momentum dependence of the self energy is weak. 
Comparing the normal and superconducting state disper- 
sion, the influence of the resonance can then be quantified 
using e.g. the relation 



15K 



V' 



50K 



-15K 



_50K 



(7) 
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FIG. 4: MDC-derived dispersions, obtained from the cen- 
tre positions of the Lorentzian fits as shown in the inset of 
Fig. [3j of the cei (red) and 02 (green) band as well as the /3i 
band (magenta) in the superconducting state at T = 15 K. 
The blue curves are the dispersions in the normal state at 
T = 50 K respectively. The insets show the real part of the 
effective self energy E'(e) = k . For temperature depen- 
dent results see Fig. 



tire Brillouin zone, the evaluation of tunnelling spec- 
tra of superconductor-insulator-normal metal (SIN) 
and superconductor-insulator-superconductor (SIS) tun- 
nelling junctions is straightforward. The SIN tunnelling 
current /(V) is given by 



/(yH£|M fc | 2 / |U(e,fc)[/(e 



eV)-f(e)}. (8) 



For incoherent (momentum non-conserving) tunnelling 
we assume a constant \Mk\ = \M \ . The results for the 
differential SIN tunnelling conductance psiN = dl/dV 
is shown in the top panel of FIG. [5] for varying coupling 
strength. With increasing coupling constant a dip-hump 
structure develops and weight is shifted into the peak 
at the occupied side. Furthermore, the hump position 
moves towards the chemical potential. Analogously to 
the cuprates, coupling predominantly influences the oc- 
cupied side [35, 



For a SIS junction and the case of incoherent single 
particle tunnelling we have with A(e) = J^k ^( e ' ^) 



I(V) = Fol- 



ds 
2^ 



A(e)A(e + eV)[f(e)-f(e + eV)} (9) 



and again calculate the differential tunnelling conduc- 
tance dl/dV, which is shown in the bottom panel of 
FIG. [5] Due to the symmetry of Eq. ([5]), the dip-hump 
feature is now strong on both sides. 

However, in both junctions we see no well pronounced 
dip for the coupling strength g 2 Xo = 1-17 x 10 3 ^|eV K 
which is the value we obtain from our comparison with 
ARPES experiments. 



It is common to plot E'(e) 



iff.. 



for a given 



reference temperature T re f well inside the normal phase. 

In FIG.2]the MDC-derived dispersions of the ot\ t 2 and 
the P\ band in the normal and superconducting state are 
shown. The real part of the effective self energies in the 
respective inset feature a broad maximum at sa 50 meV 
as well as a peak at lower energies. This peak appears 
due to coupling to the resonance at energy + files' 
where A*2 is the gap at the corresponding Fermi surfaces 
that are nested by an antifcrromagnetic wavevector. The 
shape as well as the absolute value of the real part of the 
effective self energy agree with experimental data 18. lfjj. 



Tunnelling spectra 

Another useful tool to obtain information on the den- 
sity of states (DOS) is the scanning tunnelling spec- 
troscopy (STS). Apart from the magnitude of the su- 
perconducting gap, also features due to the inter- 
action with spin excitations can be expected [l7| . 
Since we calculate the spectral function in the en- 



Discussion 

From the knowledge of the experimental value of xo 
it is possible to extract the coupling constant g. This 
value is not known for Bai_ a; K a ;Fe2As2, however it has 
been determined for optimally doped BaFei.85Corj.i5As2 
@. In order to obtain a rough estimate of the cou- 
pling constant g, we insert this experimentally deter- 
mined value xo = (3.8 ± 1.0) x lOVicV^K and ob- 
tain g w 0.15 — 0.2 eV. We also have performed extended 
calculations of the temperature dependence of the self 
energy effects discussed above, which are presented in 
the appendix. From our results we conclude that the ex- 
perimentally observed modifications of the electronic dis- 
persion when entering the superconducting state can be 
explained within a model where a continuum of strongly 
incoherent high-energy spin fluctuations provide the pair- 
ing interaction, whereas the low-energy spin-fluctuation 
modes lead to the observed electronic dispersion features. 
The coupling strength that is required to explain the ex- 
perimental data is large enough to allow for a magnetic 
pairing scenario in iron-pnictide superconductors. 
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FIG. 5: STS tunnelling spectra: Differential tun- 
nelling conductance for different coupling constants g 2 Xo = 
0.68.. .2.42 x lOVIeV K for the case of incoherent tunnelling 
at T — 15K. Units are e|Afo| 2 for an SIN junction (top) and 
e|Tb| 2 for an SIS junction (bottom). 



Methods 

We employ a tight-binding fit in orbital basis that was ob- 
tained from the DFT band structure of BaFe2 As2 by Ref. [26( |. 
The largest contribution to the density of states comes from 
the five Fe-d-orbitals in the energy region we consider, which 
allows to reduce the Hamiltonian 



H = d L(ka) [f mn (fc) + S„ 



d n (ka), 



bital m, where m — 1 ... 5 corresponds to the five orbitals 
d xz , dyz, d xy , d x 2_y2,d 3z 2_ r 2. The parameters £ m „ and e„ 
are listed in Ref. [2g|. The chemical potential at zero tem- 
perature in the undoped state is at /i c = 0. The canonical 
transformation d^(ka) = a ™(fc) dl n (kcr), where a™(fc) = 
(mk\fj.k), diagonalises the Hamiltonian, leading to eigenvalues 
£ M (fe) and eigenvectors a™(k), where /x represents the band 
index. We assume that the set of five eigenvectors for each k is 
orthonormal, Yl m a fT(fc)* a IT(fe) = &t±v, and use completeness, 
a™(k)a^(k)* = 5 mn . In order to simulate hole doping we 
apply a rigid shift of the chemical potential by an amount of 
S/j, c — — 50meV. This leads to the appearance of additional 
hole pockets around (ir,ir,kz) that favour an s ± -state [2t| . 

The unperturbed Green's function is diagonal in band in- 
dex, with normal (diagonal) and anomal (off-diagonal) com- 
ponents G^°'(e, fc) and Fji \e,k). The renormalised Green's 
functions, G M „ and F^, are not diagonal in band index due to 
interband interactions introduced by spin fluctuations. The 
Green's functions in an orbital basis, G mn and F mn , are re- 
lated to those in a band representation, G M „ and F^ v , by 

Gmn(e, k) = q™ (fc) a" (fc) * G M „ (e, k) , 
F mn (e,k) = ^a™(fe)a"(-fe)F M 4e,fc). 

Note that due to £ mn ( — k) = £ mn (fc)* the eigenvectors can 
be chosen such that q™(— fc) = q™(fe)* holds. 

Magnetic excitations can be described in terms of the 
dynamic spin susceptibility Xnj?i where the measured sus- 
ceptibility vis the sum over all orbital contributions \ — 
Y2mn XrirT [2al- Since the magnetic susceptibility is in partic- 
ular enhanced for intra-orbital coupling [28|], we assume that 
the elements Xnn of the susceptibility dominate, neglect in 
our model all inter-orbital contributions, i.e. \ ~ SnX™, 
and take only those orbitals into account where the coupling 
is most pronounced (n = xz,yz,xy). The neglected contri- 
butions have a small spectral weight in the considered energy 
range oj < lOOmeV, and are thus nearly unaffected by su- 
perconductivity. Taking them into account leads to a simple 
energy and momentum independent renormalisation of the 
low-energy electronic spectra (see below). In order to de- 
termine the relative weight between w xz — w yz and w xy we 
studied the behaviour of the electronic dispersion near the 
X point of the Brillouin zone. We find that experiments are 
best reproduced when the weight is roughly equal. We thus 
performed our calculations for w xz = w yz — w xy — wt and 
Xt z = Xt ~ Xt* ~ Xt- In order to determine the ratio 
between \t and wt we employ the total moment sum rule 
for the dynamical structure factor, as the temperature vari- 
ation of the total moment is negligible. With the above as- 
sumptions, we can separate the sum-rule into its orbital parts. 
Employing the fluctuation dissipation theorem that connects 
the dynamical spin structure factor with the susceptibility, 
S(q,w) = 2?i£ n x"c,n(<7,o;)/(l-e-'-V fc * T ), the ratio w t / X t 
is determined by the equality 



/ dw [ 

J -2A J 



d qS r , n (q,u) 



dw Id qS c ,n(q,u) 



to this basis. Here dl n (ka) creates an electron with mo- 
mentum k and spin a [we denote k — (fc, a)] in the or- 



-2A J J -2A 

with S c(r) , n (q,u) = 2fiXc(r), n (q.")/(l - e- ft " /^ T ) an d 
temperature dependent order parameter A oc y/l — T/T c . 
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We couple electrons to the spin fluctuation spectrum with 
an energy and momentum independent coupling constant g 
(instantaneous and local coupling). The retarded diagonal 
and off-diagonal self-energies are then given by 



1 y 



A <f> R 



written in terms of retarded (R) and Keldysh (K) Green's 
functions, 

v-R _ ' 2 ( n K R i n R K \ 

Z-'n — y-'nn * Xn + "nn * Xn J i 

*n = -gfl [Kn*Xn + Kn*Xn ) , 

with (A * B)(e,k) — '^2 ujg A(e — u, k — q)B(u,q) as ex- 
plained in Ref. [ljj], and with the diagonal susceptibility, 
Xn = Xnn- I n the normal and superconducting state the sus- 
ceptibility is given by \n = Xc an d Xn = Xsc respectively, 
where n = 1,2,3. We calculate the convolutions numerically 
by fast Fourier transform, using bare Green's functions (a pro- 
cedure supported by the numerical studies in Ref. [33) with 
a broadening parameter 8 = 4meV. 

Using the completeness relation = 5 mn , we 

can write the Dyson equation for the renormalised retarded 
Green's function in Nambu-Gor'kov space, 

G mn (e, fc) = G^nn E mn = ^ ^ ®^ ^ ^ ? 

in the following way 

Gmn (6) k) = ^ a"(fc)q^(fc)* 

x [(Z n (e, k) (e + - ^(e, fc)a z - A n (e, fc)*,] . 

where 1 and &i are the 2x2 unit matrix and the Pauli matrices 
in Nambu space, respectively. The renormalised dispersion 
and order parameter as well as the renormalisation function 
are given in terms of the retarded diagonal and off-diagonal 
self energies, 



£(e,k) = C M (fc) + 



E*(e,fc) + E*(-e,-fc)* 



A n (e,fc) = A fe + <f>*(e,fc), 

, £«(e,fc)-E*(-e,-fc)' 



Z n {t,k) 



2(e + io) 



using the fact that 4>^(e, k) = $^(— e, — fc)*. 
The equation for G^n ( e , fc) is inverted numerically to obtain 
G^ n (e, fc). The spectral function is then obtained from the 
retarded Green's function via 



For the numerical calculations we use a 2 x 2 x 2 k-mesh 
and 2 7 points in energy space. 

A high-energy cutoff of lo c — 200 meV was introduced in 
the spectrum of spin excitations. The exact value of this cut- 
off, as well as spectral weight beyond this cutoff and the pre- 
cise variation of the spin susceptibility in the vicinity of the 
cutoff energy are however not of importance, as any change 
in the high-energy contributions only leads to an additional 



renormalisation factor Z HE for the low-energy electronic ex- 
citations. Thus as the high-energy part of the susceptibility 
varies with temperature, so does Z HE . Taking the dispersion 
relation Eh and the imaginary part of the effective self energy 
E"(fe) as obtained from Eq. ((6| in the main text, the spectral 
function in the superconducting state can be well approxi- 
mated by 



AJk) = --1m 

TV 



1 



Z € {k)(e + i5) -e fc + 



A g (fc) 2 



where A e (fc) is the renormalised order parameter. The 
renormalised value of the gap is obtained by considering 
the quantity A e (k)/Z e (k). The total renormalisation fac- 
tor Z t (k) can be expressed as a product of high- and low- 
energy renormalisations, Z e (k) — Z HE (k)Z LE (fc), with an 
energy-independent (on the low energy scale) Z HE (k), and 
with Z LE (k) = 1 - E"(fe)/(e + i5). In order to compare the 
results for different temperatures we have to include the high 
energy contributions to the renormalisation of the electronic 
dispersions, i.e. e h >em = Eh/Z HE (fc). The high-energy renor- 
malisation factor is only weakly momentum dependent due 
to the weak momentum dependence of the susceptibility at 
high energies, and consequently we neglect this momentum 
dependence in our calculation. We are able to reproduce the 
experimentally observed high-energy part of the dispersions 
with Z HE {T = 50K) = 1, and determine Z HE for the other 
temperatures so that the dispersions merge the normal state 
dispersion at T — 50K for high energies. For the tempera- 
tures T = 15 — 35K Z HE varies slowly between 0.9 — 1 as 
it is shown in Fig. [S] Its value is < 1 to account for the 
overestimate of x" near uj c at lower temperatures. 



0.6 



20 



25 
T[K] 



30 



35 



FIG. 6: Temperature dependence of the high energy renor- 
malisation factor Z HE corresponding to a high-energy cutoff 
for the spin excitation spectrum of u c = 200 meV. 
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TEMPERATURE DEPENDENCE OF THE SPIN FLUCTUATION SPECTRUM AND THE RES- 
ONANCE MODE 

We model the spin fluctuation spectrum by the procedure explained in the main article. The functional form of 
x"(w,q) is given by Eqs. (1) and (2) in the paper. In FIG. [7] the temperature dependence of the bosonic mode is 
shown. We present the energy dependence of the dynamic spin susceptibility at fixed momentum q = Q, where Q 
is the antiferromagnetic wave vector, as well as the energy-momentum dependence of the spectrum. With decreasing 
temperature, spectral weight of the two-particle excitation spectrum is shifted into the energy region |e| < 2A(T), 
leading to a well pronounced resonance feature that is sharp both in momentum and energy. Near T c the resonance 
disappears and the energy dependence obeys the same shape as the normal state dispersion at T = 50 A' . 



T=13.1K 



T=31.0K 






0.02 0.04 I) .06 




T=18.2K 



T=33.6K 







T=23.6K 



T=35.5K 






0.02 0.04 0.06 




T=27.7K 



T=36.6K 







FIG. 7: Dependence of the spin fluctuation spectrum on temperature. On the left of each panel the shape of the 
dynamic spin susceptibility at the antiferromagnetic wave vector, x'(ui,Q), is shown (red line). The blue curve corresponds 
to the normal state spectral function at T — 50K. On the right hand side the energy and momentum dependence around the 
wave vector Q is presented. 
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INFLUENCE OF THE BOSONIC SPECTRAL SHAPE ON THE ELECTRONIC SPECTRA 

In addition to the results in the paper, which are given for the temperatures T = IbK and T = 50K, we performed 
the same calculations for various temperatures (T = 6.8 — 36.6K). The appearance of the bosonic resonance leads 
to an effect on the electronic dispersion which is characterised by the development of a peak in the real part of the 
effective self energy as well as a hump feature in the imaginary part of the effective self energy. Both effects clearly 
show the same temperature dependence as the resonance in the bosonic spectrum as presented in FIG [8] and FIG [H 
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FIG. 8: Temperature dependence of the real part of the effective self energy at the /3i-band. With increasing 
temperature, the self energy shows the development of a broad hump and a peak appearing at low energies. 
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FIG. 9: Temperature dependence of the imaginary part of the effective self energy at the ai-band. We present the 
results for each temperature (red line) compared to the normal state imaginary part of the effective self energy at T = 50K 
(blue line). The black curve presents the difference of both, scaled by a factor of 3. 



